Raman modes of the deformed single-wall carbon nanotubes 
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With the empirical bond polarizability model, the nonresonant Raman spectra of the chiral and 
achiral single-wall carbon nanotubes (SWCNTs) under uniaxial and torsional strains have been 
systematically studied by ab initio method. It is found that both the frequencies and the intensities 
of the low-frequency Raman active modes almost do not change in the deformed nanotubes, while 
their high-frequency part shifts obviously. Especially, the high-frequency part shifts linearly with 
the uniaxial tensile strain, and two kinds of different shift slopes are found for any kind of SWCNTs. 
More interestingly, new Raman peaks are found in the nonresonant Raman spectra under torsional 
strain, which are explained by a) the symmetry breaking and b) the effect of bond rotation and the 
anisotropy of the polarizability induced by bond stretching. 

PACS numbers: 63.22.-|-m, 78.30.Na 
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I. INTRODUCTION 

Single-wall carbon nanotube (SWCNT) is a kind of 
nanoscale molecule obtained by wrapping a graphene 
sheet into a seamless cylinder. Since its discovery- in 
1991 by lijima, the SWCNT has been of great interest 
in mesoscopic physics and nanotechnology because of its 
exceptional electronic and mechanical properties^. 

At the same time, much attention has been 
devoted to the vibrational properties of the 
SWCNTsM^>^2^^ao^i2^a4^a6^ The excita- 

tions from the electronic density-of-state (DOS) peaks 
below the Fermi level to those above it are enhanced 
greatly due to the existence of the Van Hove singularities 
(VHS)^ caused by the special one-dimensional (ID) 
structure of the SWCNTs. As a result, the Raman 
spectroscopy has been considered as a powerful tool 
to detect the diameter-selective phonon modes of 
SWCNTsi>i^>ia>20. The low-frequency peaks in the 
Raman spectrum are attributed to the radial breathing 
modes (REM, R band), where all carbon atoms are 
subject to an in-phase radial displacement. It was found 
that the REM frequency depends only on the tube 
diameter in an inverse proportion to it, but not on its 
chirality. On the other hand, the higher frequency peaks 
in the Raman spectrum are caused by the out-of-phase 
vibrations of the neighboring carbon atoms parallel to 
the tube surface (tangential modes, T band), which are 
related to the i?2g(2) phonon mode at ~1580 cni^^ in 
graphite. In the achiral tubes, one can distinguish two 
tangential modes with atomic displacements parallel and 
perpendicular to the tube axis, respectively. 

Eeyond traditional force-constant method^, the state- 
of-the-art first-principle methods have already been 
used to calculate more accurately phonon dispersion 
curves^i^i22ii^'2^, which offers a unique technique to in- 
vestigate the vibrational properties of the nanotubes or 
nanoropes. Moreover, some of the first-principle meth- 
ods can be extended to predict nonresonant Raman 



intensities2£j which, however, needs very intensive cal- 
culations and is less accurate than the purely vibrational 
calculations, leading to proposal of a simple empirical 
bond polarizability model to calculate the Raman inten- 
sities of the fuUerene o^^'^^ and nanotubes*. It is proved 
that the method can provide reasonable nonresonant Ra- 
man spectra in the most cases, including the nonresonant 
Raman spectra of the finite- length SWCNTs^^. 

Recently, the mechanical deformation effect on the 
SWCNT electronic properties has been successfully stud- 
ied by using a simple picture based on the tt elec- 
tron tight-binding approximation^. Some interesting 
behaviors are found, e.g., a symmetry breaking due 
to the deformation may remove the degeneracy of en- 
ergy levels. Experimentally, resistance of SWCNT tran- 
sistors was found to vary significantly under bending 
and stretching''^, and an observed torsion of a metal- 
lic SWCNT was speculated to open a small band gap 
of the tube^. Several theoretical studies on the uniaxi- 
ally and torsionally deformed SWCNTs^OiSS^^^ Yi&Ye 
shown that their band structures are determined by their 
chiral symmetries and the kinds of strains. Most of the 
deformed SWCNTs show a metal-semiconductor transi- 
tion, occurring repeatedly with increasing strain. It is 
known that under uniaxial strain, the derivative of band 
gap over strain is the largest for zigzag tubes and de- 
creases with increase of the chiral angle. In contrast, un- 
der torsional strain, the derivative becomes the largest 
for armchair tubes and decreases as the chiral angle de- 
creases. In particular, armchair tubes under uniaxial 
strain and metallic zigzag tubes under torsional strain 
remain to be metallio^. 

At the same time, the deformation effect on 
the SWCNT vibration properties, especially its Ra- 
man spectroscopy, has attracted much attention in 
expcrimcnt a''''i''*'^^''*°i'^^'^^i'^^ and theoretical studiesi^i^. 
Kahn et alr^ studied the hydrostatic pressure depen- 
dence of the Raman modes for different nanotubes, and 
found that their shifts under pressure differ from the 
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corresponding: ones of graphite. It is found in Refs. 
I37ll4^l43l that the RBM shifts towards higher wavenum- 
bers with increasing hydrostatic pressure at a rate of 
7~10 cm~^GPa~^. With a generaUzed tight-binding 
molecular dynamics (GTBMD) scheme, Venkateswaran 
et alJ^^ analyzed the effect of the intertube Van der Waals 
interaction on the R band shift. And Refs. [s^ andlH re- 
ported that before a structural transition at near 2 GPa, 
the tangential modes exhibited essentially the same pres- 
sure dependence of 5.7-5.8 and 5.3-6.1 cm~^GPa~^, re- 
spectively. More details about this hydrostatic pressure 
effect on the Raman modes can be found in the review 
article^. Nevertheless, all of these researches focuses on 
the hydrostatic pressure effect, and the uniaxial and tor- 
sional strain effects have not been studied. 

In this work we will study systematically the nonreso- 
nant Raman spectra of chiral and achiral SWCNTs un- 
der uniaxial and torsional strains. The Raman intensity 
is calculated by using the eigenvectors of the vibrational 
modes obtained by ab initio method, and the polarizabil- 
ity parameters are obtained by using the empirical bond 
polarizability model^i^^. The detail of methodology is 
presented in Sec. II. The numerical results are discussed 
in Sec. III. Finally, this work ends with some conclusions 
in Sec. IV. 



II. MODEL AND METHOD 

A. Cumulant force-constant method 

In most of published papers on the phonon disper- 
sions of the SWCNTs obtained by ah initio calculations, 
the low-frequency dispersions were not well reproduced 
mainly due to use of the interatomic force-constant (IFC) 
method. In Ref. [H, by means of Cumulant Force Con- 
stant (CFC) method, L. H. Ye et al. successfully calcu- 
lated phonon dispersions of the nanotubes and graphite. 
Details about comparison between the CFC and IFC 



methods are available in their original paper—. In this 
work, we use CFC method to construct the dynamical 
matrices for both of the undeformed and deformed SWC- 
NTs, and then calculate their phonon dispersion curves. 

First, we summarize the reason of employing CFC 
method in this work: 

a) . The ID-like geometrical structure of the car- 
bon nanotube makes its physical properties highly 
anisotropic, for which it is definitely necessary to use the 
CFC instead of the IFC. Comparatively, the CFC is also 
an approximated method, but it keeps the most rota- 
tional symmetries, which means although the numerical 
values of the dynamical matrices are obtained approx- 
imately, the symmetries of the dynamical matrices are 
correct. This is the reason why it works so well for an 
anisotropic system. 

b) . In the IFC method, in order to minimize the trun- 
cation effect, one should use a large supercell, spending 
much computational time. At the same time, the nu- 
merical noise in force calculations may become very se- 
rious if too large supercell is used, which may reduce 
the accuracy of phonon calculation. On the other hand, 
because the CFC method does not need to perform an 
explicit truncation on the interatomic interactions, as 
shown in Ref. [2^, so, only a very small supercell (1x1x2) 
in the CFC method is sufhcient to calculate accurately 
the phonon modes on the F point and the X point at Bril- 
louin zone boundary. This also means that the phonon 
modes of even achiral SWCNTs, such as (12, 4) tube con- 
taining 208 atoms in its one unit cell, can be efficiently 
calculated under current hardware conditions. 



B. Empirical bond polarizability model 

The time-averaged power flux of the Raman-scattered 
light at a given direction of a solid angle (Kl in a fre- 
quency range between uj f and ujf + dojf is related to the 
differential scattering cross section as follows. 



ijj}uJi[{n{uj)) +l]h-K VaVpHa'yl3\{l^)w^W\, (1) 
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(j) = ^^(5n|j). (3) 

no 



, , v-^ * / -N /-N 1 rr/ \ r/ M Here, uJi and wr are the frequencies of the incident 

"TPA K , a7 KJ) PA yJi 2(^-^ " ^ ■'^J'and scattered light; w = - w/ is the Raman shift, v 

(^2) and w are the corresponding polarization unit vectors of 
the incident and scattered light, respectively, {n {uj)) is 
with the Bose factor, M„ is the mass of the nth atom, and 
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iVj and {5n \ j) are the frequency and {Sn) component of 
the jth mode. The coefficient tt^^ ^ in Eq. ^ connects 
the polarization changes to the atomic motions^, which 
is obtained by expanding the polarizabihty tensor tt" in 
terms of the atom displacements u^, with 



Adopting the nonresonant empirical bond-polarization 
theorySi2Ij£§.^ the Eq. @ can be rewritten as: 
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{Sa^yT-p +6(3jra-2rafpf^). (5) 



Here, f is the unit vector of the vector r connecting the 
n and m atoms linked by a bond, ay and a± represent 
the static longitudinal and perpendicular bond polariz- 
abihty, respectively, and aj| = (^^^ , a'± = (^)o- 
The values of , a_i_, a'| and a'j_ are given empirically as 
a function of the bond length between two carbon atoms. 
For undeformed SWCNTs, there is only single type of 
bonds and the bond polarizabihty model is completely 
defined by three parameters^: a — 2a'| -I- a'^ = 4.7 A^, 

P = a\^-a'j_= 4.0 A2, and 7 = a|| - a_L = 0.04 A^. It 
is known, however, that the polarizabihty parameters of 
carbon are similar for a variety of carbon materials. Fur- 
thermore, the relative intensities of the Raman modes are 
not so sensitive to small changes of the bond polarization 
parameter values, except for the lowest E2g mode. Thus 
in the present work, we still use the above polarization 
parameters for the deformed SWCNTs, which can only 
lead to a small error in the peak height, and does not af- 
fect the peak positions. In fact, Guha et aL— ^ chose differ- 
ent bond polarizabihty parameters to calculate the Ra- 
man spectra of the fuUerenes, causing only small changes 
of the relative Raman intensities, which is also found in 
our calculations. For example, even use of a set of pa- 
rameters obtained by interpolating those of the fuUerenes 
and nanotubes does not bring any qualitative changes of 
our results. 

In all our calculations, the line shape of each peak is 
assumed to be Lorentzian and its line-width is fixed at 
1.68 cm^^. The tube axis is taken along the Z direction, 
and the X and Y directions are taken to be perpendicu- 
lar to it. In Ref. 8^ two possible geometries for the light 
polarization are mentioned: the VV and VH configu- 
rations. In the VV configuration, the incident and the 
scattered polarizations are parallel to each other, while 
they are perpendicular to each other in the VH configu- 
ration. After careful calculations, we can find all of the 
Raman-active peaks in the VV configuration. For the 
VH configuration, there is no any qualitative change of 
our results compared with those obtained from VV con- 



figuration. Due to this reason, and also because we pay 
our main attention to the shift of the Raman mode fre- 
quency under applied strains, only the results in the VV 
configuration have been presented in this work. 



C. Calculation details 

Our first principle calculations are based on the total 
energy plane-wave pseudopotential method in the local 
density approximation (LDA) . The Vienna ab initio sim- 
ulation package (VASP)— 1^ is employed in this paper. 
In the most recent version of VASP, the interaction be- 
tween the ions and electrons is described by the projec- 
tor augmented-wave^" (PAW) in the implementation of 
Kress and Joubert^^. The PAW method is always con- 
sidered to be able to give a more accurate and reliable 
result than the ultrasoft potential. In the present PAW 
potential, the 2s and 2p orbitals are treated as valence 
orbitals, and a plane- wave cutoff of 400 eV is used to 
obtain reliable results. 

A supercell geometry was adopted so that the nan- 
otubes are aligned in a hexagonal array with the clos- 
est distance between adjacent nanotube surfaces being 
at least 10 A, which is found to be large enough to pre- 
vent the tube-tube interactions. 

The K-points sampling in the reciprocal space is a uni- 
form grid (Ixlxn) along the nanotube axis {Z direc- 
tion) in our calculations, the maximum spacing between 
k points is 0.03 A~^ and the Gaussian smearing width 
is 0.03 eV. This means that n is obtained by dividing 30 
with the supercell length along the tube axis. For exam- 
ple, for a supercell (1x1x2) constructed by the unit cell 
of a zigzag tube, the n equals to 30/ {T^ig * 2) w 6, where 
Tzig ~ 2.451 A is the unit cell length of the zigzag tube. 
The increasing K-points sampling has no significant ef- 
fect on the frequencies of the final Raman active modes, 
e.g., doubling the number of special k points only causes 
less than 4 cm~^ change on the most F point phonon 
frequencies. 
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All SWCNT geometrical structures used in this work 
have been fully optimized before performing further force 
constants (FCs) calculation. The structural optimization 
process is almost the same as Ref. [1^ Finally, the opti- 
mal structure is obtained when the residue forces acting 
on all the atoms were less than 0.005 eV/A for achiral 
SWCNTs and 0.01 eV/A for chiral SWCNTs. The equi- 
librium lattice constants and averaged nanotube radii are 
listed in Table I. Our structure parameters of (5, 5), (10, 
0) and (10, 5) tubes are almost the same as those in Ref. 
[Sa . which are also obtained by VASP. 

Due to the high symmetry of an ideal SWCNT, so there 
exists only one inequivalent atom. Thus the FCs of one 
atom are sufficient to construct the force constant matrix 
of the whole tube by applying the tube symmetry. 



TABLE I: Optimized structural parameters of the unde- 
formed nanotubes. -Ro is the average radius, and To is the 
lattice constant along the tube axis. 



Tubes 


Ro{A) 




(6, 0) 


2.397 


4.230 


(7, 0) 


2.773 


4.235 


(8, 0) 


3.159 


4.233 


(9, 0) 


3.541 


4.237 


(10, 0) 


3.927 


4.237 


(11, 0) 


4.313 


4.238 


(5, 5) 


3.408 


2.448 


(6, 6) 


4.075 


2.448 


(10, 5) 


5.173 


11.21 


(12, 4) 


5.636 


15.28 


III. 


NUMERICAL RESULTS AND 
DISCUSSIONS. 





In this section, the Raman spectra of the infinite-length 
deformed SWCNTs are investigated. Firstly, the tensile 
strain effect will be discussed in Sec. Ill A, and the tor- 
sional strain one will be given in Sec. Ill B. 



A. Tensile strain 

The phonon dispersion relations and density of states 
for zigzag (10, 0) tube are shown in Fig. 1, from which, 
one can find that the low-frequency dispersions are well 
produced, and when the tensile strain is applied, the low- 
frequency vibration modes almost do not change, while 
the high-frequency vibrations show sensitive response to 
the tensile strain. 

These characters can be seen more clearly in the Ra- 
man spectra. First we will consider the Raman spectra of 
zigzag tubes as an example. Because the cylindric struc- 
ture and the neighboring relationships do not change un- 
der the tensile strain, we find the intensities of the Ra- 
man active modes almost do not change at all. Thus only 



the frequency shifts for the Raman active modes under 
tensile strain are given in Fig. 2. In the present work, 
the tensile strain is defined as e = (T/Tq) — 1, where 
T and Tq are the lengths of the axial unit cells of the 
deformed and undeformed nanotubes, respectively. The 
values of Tq have been given in Table I. After applying 
tensile strain, only the positions of the carbon atoms are 
relaxed to obtain a reliable optimized structure. It can be 
seen clearly from Fig. 2 that when the tensile strength 
increases: a) The frequencies of the low-frequency Ra- 
man active modes almost do not change. While the high- 
frequency part of the Raman spectra is shifted downward, 
which consists with the available experimental and theo- 
retical resultsi^>27i38^do.4i^i43^. b) The deformation 
does not change the number of the Raman active vibra- 
tions, c) The shifting rate of the Raman modes is almost 
the same for the nanotubes with different radii, i.e., they 
shift linearly. 

In fact, above three results are reasonable and can be 
understood as follows. Because different nanotubes have 
almost the same Poisson ratio a of 0.2, the tube radius 
only changes about 1%, even if the tube is elongated by 
5%. And for lower Raman modes, the vibrational fre- 
quencies are mainly affected by the tube diameter, mak- 
ing their peak positions do not move much. Especially, 
for the radial breathing mode (RBM), which is found 
to be inversely proportional to the tube diameter, only 
very small shift (about 1 cm~^) is found when the tensile 
strength increases up to 5%. 

But the bonds parallel to the tube axis will be elon- 
gated when tensile strain is applied. For higher Raman 
frequencies, the vibrations are sensitive to the local bond 
structure, especially its length. When bond length is 
elongated, the interaction among nearby atoms will be- 
come weaker, making the Raman frequency smaller. 

Now we would like to discuss why the number of Ra- 
man active modes does not change under the tensile 
strain. Following the discussion of Ref. 52, the num- 
ber of Raman active modes can be determined by using 
the nonsymmorphic rod-groups of SWCNTs. All achiral 
carbon nanotubes can be shown by this method to pos- 
sess only 8 Raman-active phonon modest, and for all 
chiral carbon nanotubes, 14 Raman-active modes can be 
found. Furthermore, for the first-order Raman scatter- 
ing, the point group is sufficient because this process does 
not change the wave vector k. Another fact is that the 
point groups of achiral (zigzag or armchair) and chiral 
SWCNTs are D2nh and Z?jv, respectively, where n is the 
tube index and N is the number of graphene hexagons 
in the nanotube unit cell. These point groups do not 
change under tensile strain. Correspondingly, simply fol- 
lowing the discussion in Ref. [s^ . one can conclude that 
when tensile strength increases, the Raman active vibra- 
tion number of zigzag nanotubes will keep 8, which has 
been shown by our numerical results. 

Then let us analyze the detail of frequency shift of 
Raman active modes. From Fig. 2, one can find that 
the modes below 1000 cm~^ is insensitive to the tensile 
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[(states/C atom) 



(a). 0% tensile 



(b), 2% tensile 



(c) 4% tensile 



FIG. 1: Calculated phonon dispersion curves and density of states are shown in the left and right panels, respectively, in each 
sub-figure. For (10, 0) zigzag nanotube under the tensile strains of 0%, 2% and 4%. 




Tensile Strain (%) Tensile Strain (%) Tensile Strain (%) 

(a). (6,0) (b). (8,0) (c). (10,0) 



FIG. 2: Frequencies of Raman active modes for (n, 0) zigzag tubes (n=6, 8 and 10) under tensile strain. 



strain. And above 1000 cm~^, the frequencies of Ra- 
man active modes decrease linearly with tensile strain, 
for which there exist two different decreasing slopes for: 
(i). Eig, and (ii). and Aig modes. 

Similar phenomenon has been found in the graphite 
system^. Following Ref. 113, we also calculated the 
Griineisen tensor elements (GTEs) k = ~ ^in(i"e) ' where 
Lo is the frequency of Raman active modes, and e is the 



tensile strain. We plot the GTEs vs the tube index n 
of the zigzag tubs in Fig. 3, and further fit them with 
a same fitting function in the same plot. The fitting 
function is obtained based on the following factors: a) 
The GTEs must have the same dependence on the tube 
radii as the frequencies of Raman modes. In Ref. [ssl . 
the higher Raman frequencies depend on the tube ra- 
dius via the confined wave vector q around the cqua- 
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FIG. 3: (Color online) The GTEs for the highest three Raman 
active modes of zigzag tubes (n, 0), (n — 6, ...,11). The 
lines are fitting results (see text), and other markers (circles, 
diamonds and squares) are GTEs obtained by first-principle 
calculations. 



torial direction, where q oc 1/r oc 1/n. Here we fol- 
lowing this idea and set the nonoscillating part of the 
fitting function as A + ^S;-, where A and B are unde- 
termined coefficients, whose quadratic form fulfils the 
phonon dispersion near the T point of graphite, b) The 
above k form still can not describe the oscillating behav- 
ior of the GTEs. Thus an extra oscillating term must be 
added into the fitting function, which can be written as 
^sin(^n -I- (jjo), where the period has been set as 3 to 
match the oscillating period of the energy gap of zigzag 
tubes, and the denominator is used to let this term 
quickly vanish when n — > oo. So the final fitting func- 
tion can be written as: k — A + ^ + ;^sin(^n -I- 0o)- 
After making a simple least-squares fit, finally we can 
obtain KE^g {r oo) = 0.897, KAig (r — > oo) = 2.05, and 
KE2 {f ~^ 00) = 1.95, which can be checked by future 
experiments. This also offers a powerful tool to deter- 
mine the strength of the tensile strain, i.e., after detect- 
ing the Raman-active modes under undeformed and de- 
formed conditions, one can easily estimate how big the 
applied tensile strain is. 

Next, let us consider the armchair and chiral SWC- 
NTs. We find those three key characters found for the 
Raman spectra of zigzag tubes are also hold for these 



SWCNTs. The reason should be the same as that for the 
zigzag tubes. Because we mainly concern the frequency 
shift, so the Raman spectra will not be shown, and in- 
stead only the frequencies of the Raman active modes for 
some armchair and chiral tubes under tensile strains are 
shown in Fig. 4, which can already represent all the key 
characters clearly. 

The GTEs for these tubes are also calculated. It is 
interesting to note that the GTEs of the high-frequency 
Raman modes of armchair tubes can also be divided into 
two groups, and their values are very close to those of the 
zigzag tubes. For example, for (6, 6) tube, KEig = 2.35, 
KA,g = 0.861, and ke^^ = 0.794. And for (12, 4) tube, 
the GTEs of the six highest frequency Raman modes are: 
2.25, 2.27, 2.28, 0.994, 0.907 and 0.846. So, we suppose 
that for any SWCNTs, there are only two kinds of GTEs 
for the high-frequency Raman modes. In fact, detailed 
analysis shows that those modes with the larger GTEs 
are the tangential with atomic displacements parallel to 
the tube axis, and those modes with the smaller GTEs 
are the tangential modes with atomic displacements per- 
pendicular to the tube axis. This fact could be clearly 
seen from Fig. 5, showing the atomic displacements of 
three high-frequency Raman-active tangential modes of 
(6, 6) tube. 

Finally, one may notice that it is different for the ten- 
sile strain to affect the electronic and vibrational prop- 
erties of SWCNTs. As mentioned in the introduction 
part, most of the deformed SWCNTs show a metal- 
semiconductor transition, occurring repeatedly with in- 
creasing strain. Correspondingly, the linear optical re- 
sponses of the deformed carbon nanotubes also show pe- 
riodical change with increasing strain^^ . But as shown in 
our calculations, the frequencies of Raman active modes 
change monotonously with applied tensile strain and 
their intensities only change very little, which is consis- 
tent with the available experiments^li^^i^daiilii^ii^iS. 



B. Torsional strain 

In this part, the torsional strain effect on the Raman 
modes is investigated. Here, three tubes with different 
chiralities are taken as an example, i.e., armchair (5, 5), 
zigzag (10, 0) and chiral (12, 4) tubes. The torsion angle 
7 is defined in a similar way to that in Ref. [30l i.e., it 
makes any point (rc, rj) on the 2D projection of the tube 
surface move to (rc — tan (7) rt,rt) (see Fig. 6). The fol- 
lowing method is used to find a new periodical structure 
under torsional strain: firstly, an origin of the coordi- 
nates is selected as the site O, from which, we can define 
a chiral vector C and the original translational vector 
Tq. Then, we find an equivalent site B for the site O. By 
projecting the vector OB to the Tq, a new translational 
vector OB' can be produced as a tube period, and the 
angle between OB and OB' is the corresponding torsion 
angle 7. After structure relaxation on both the lattice 
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Tensile Strain (%) Tensiie Strain (%) Tensile Strain {%) 

(a). (5, 5) (b).(6,6) (c).(12,4) 



FIG. 4: Frequencies of Raman active modes for some armchair and chiral SWCNTs under tensile strain, (a) (5, 5), (b) (6, 6), 
(c) (12, 4). 



o)„ =1407.4361 cm"' 



CD„ =1535.3345 cm"' 



cBj =1647.7134 cm"' 



FIG. 5: (Color online) The atomic displacements of three 
Raman-active tangential modes of (6, 6) tube. tJo is the fre- 
quency of Raman-active modes of undeformed tube. 



constant along the tube axis and the atomic positions, 
we can finally obtain a reasonable structure which can 
be used for further first principle calculations. And in all 
of our calculations, the system keeps a cylinder structure 
and does not appear any structure transformations^^ or 
defects^. 

The Raman spectra of different SWCNTs under tor- 
sional strain have been shown in Fig. 7. Because of the 
symmetry, we only consider the case of positive 7 for the 
armchair and zigzag tubes. Compared with the situation 
under the tensile strain, we can find two similar key char- 
acters: a) Not only the frequencies of the low-frequency 
Raman active modes, but also their intensities almost do 
not change. Some of the high-frequency modes in the 
Raman spectra are red-shifted and some ones are blue- 
shifted, b) The Raman active vibration number changes 
under torsional strain. 




FIG. 6: (Color online) 2D projection of the nanotubes under 
the torsional strain, where 7 is the torsion angle and C is the 
chiral vector. To and Tnen, are the translational vectors for 
the undeformed and the deformed tube. 



In fact, the first character can be explained by the same 
reason as that under tensile strain. For example, the radii 
of the optimized (12, 4) tubes changes only about 0.2% 
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FIG. 7: The Raman spectra of different SWCNTs under torsional strain, (a) (5, 5), (b) (10, 0), (c) (12, 4). 



when the torsional stain is increased up to 6 ° . And the 
two bonds mirroring each other with symmetry (see 
Fig. 1 of Ref. i52i) will be elongated and shortened un- 
der the torsional strain, respectively, making thus some 
Raman modes red-shifted, and some others blue-shifted. 
The second character may be not so obvious like the first 
one. In fact, this character manifests itself in two dif- 
ferent ways, originating from different mechanisms. The 
frequencies of Raman active modes for the three tubes 
under torsional strains are shown in Fig. 8, from which, 
it is seen that in the range from 1400 to 1600 cm"-'^ for 
all the tubes, some modes split into two new modes, ob- 
viously arising from symmetry breaking. Since the tor- 
sional strains applied in this work is smaller, the two new 
modes are very close to each other. 

On the other hand, some new low-frequency modes 
(about 400 cm-i) appear in the zigzag tubes, which can 
not simply explained by symmetry breaking. In order to 
find the origin of these new modes, the separate contri- 
butions to the Raman intensities from the three different 
terms in Eq. ([5]) are investigated. The same technique 
had been used in Ref. [H to gain a deeper insight into the 
significance of the separate contributions due to different 
hydrocarbon parameters and the fit parameters for Cgo- 
Only the results under torsional strain are given in Fig. 9 
because the three parts vanish in undeformed SWCNTs. 
From Fig. 9, one can find that the torsional strain do 
break the forbidden rule of these modes^ - and affect their 
Raman intensity dramatically. With torsional strain, the 
contributions from the second and the third term in Eq. 
([5]) enhance greatly. The second term in Eq. (O corre- 



sponds to the contribution from the anisotropic part of 
the polarizability under bond stretching, and the third 
term corresponds to the contribution from bond rotation. 
Obviously, the three bonds are elongated differently un- 
der torsional strain, which means the anisotropic effect 
will be important. And because torsional strain makes 
the bonds rotate, the third term will also enhance. Al- 
together, the scattering cross section and thus the total 
Raman intensity at a specific frequency will increase un- 
der torsional strain. 



Finally, the GTEs of high-frequency Raman modes of 
these tubes under torsional strain are calculated, which 



are defined in the present work as Hf 



16 



i In (J 



3 d taii'^ 7 



Here, 



uj is the frequency of Raman-active modes, and 7 is the 
torsional strain. This definition is obtained by assum- 

r -1 3Kt 

ing = which is taken from Ref. \54. The Kt 

for the 4 highest Raman- active modes of (5, 5) tube are 
found to be -26.73, -26.52, 19.85 and 15.65, respectively, 
and Kt for the 5 highest Raman-active modes of (10, 0) 
tube are -13.68, -30.68, -28.37, 17.95 and 19.80, respec- 
tively. Because of the complexity of the high-frequency 
Raman modes of (12, 4) tube, here only its two most in- 
tense modes at about 1600 cm^^ are considered, whose 
Kt are -32.29 and 14.03. One can find that these values 
do not obviously relate to each other, which is different 
from those under the tensile strain. So, more investiga- 
tions should be done to find the intrinsic mechanism to 
cause the shifts of the Raman-active modes under tor- 
sional strain. 



9 



200<>- 



3 — O 



3 — O 



-e^-f^^-o — o - 
o o o o o — o 



1 2 3 4 5 

Torsion Strain (Degree) 
(a). (5, 5) 



g 800 



8 § 6 



coo — e e- 



o o o — e e- 

o o ~e — s e- 



o o o — e e- 



12 3 4 

Torsion Strain (Degree) 
(b).(lO. 0) 



5 800 
I 

600 



O 



-e — e Goo 



-e — e 



-o — o o o o 



8 8 8 



o o 



-2 2 

Torsion Strain (Degree) 
(c).(12. 4) 



FIG. 8: Frequencies of Raman active modes vs torsional strains for three SWCNTs. (a) (5, 5), (b) (10, 0), (c) (12, 4). 
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FIG. 9: (Color online) Raman peak intensities of deformed 
(10, 0) tube with torsional angle of 2.2047 ° , contributed from 
the three individual terms in Eq. ([5}. The asterisk and pen- 
tagram indicate the positions of new Raman-active modes. 
Panel (I), (II) and (III) give the contributions from the first, 
the second and the last term in Eq. ((SJ, respectively. Notice 
that the total intensity is not the sum of those in the three 
panels due to interference effects between them. 

IV. CONCLUSIONS 

In this paper, the nonresonant Raman spectra of chiral 
and achiral SWCNTs under uniaxial and torsional strains 
have been systematically studied by the ab initio calcu- 
lations and the empirical bond polarizability model. The 



cumulant force-constant method is used to construct the 
dynamical matrices. Under tensile strain, it is found that 
the frequencies and the intensities of the low-frequency 
Raman active modes almost do not change, while their 
high-frequency part shifts downward linearly with two 
different slopes for any kind of SWCNTs. And the Ra- 
man active vibration number does not change under the 
tensile strain. Under the torsional strain, the frequen- 
cies and the intensities of the low-frequency Raman ac- 
tive modes still change very little, but some of the high- 
frequency modes are red-shifted and some others blue- 
shifted. More importantly, new Raman peaks are found 
in the nonresonant Raman spectra under torsional strain, 
which are explained by a) the symmetry breaking and b) 
the effect of bond rotation and the anisotropy of the po- 
larizability induced by bond stretching. 
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